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Abstract 

The prospects for a precise and model independent determination of \ V u b\ from 
inclusive and exclusive semileptonic B decays are reviewed. 



The next generation of B decay experiments will test the flavor sector of the standard model 



<N 1 I. INTRODUCTION 

> 

o 

CO 

\D ! with hieh precision. The basic approach is to determine the elements of the CKM matrix using 
different methods and then check for the consistency of these results. In practice this amounts 
to determining the sides and angles of the unitarity triangle from CP conserving decays and 
from CP asymmetries, respectively. For these checks to be meaningful, a precise and model 
Qh! independent determination of the b — > u CKM matrix element, IV^I, is very important. 

CP violation has only been observed in kaon decay arising from K° — K° mixing. This 
observation can be accommodated in the three generation standard model using the otherwise 
free parameter 5 in the CKM matrix, but this description of CP violation has not yet been 
tested. Moreover, to explain the baryon asymmetry of the universe, other sources of CP 
violation are needed [|lj . Many extensions of the standard model have new particles with weak 
scale masses, and could give large contributions to CP asymmetries or flavor changing neutral 
current processes (like K° — K° mixing, B° — B° mixing, B — > K*^, etc.). 

At the present time, the magnitude of V u b is determined by comparing experimental results 
on the endpoint region of the electron spectrum in inclusive B decays with phenomenological 
models , or by comparing experimental results on B — > p £ v and B — » 7r £u with phenomeno- 
logical models and lattice QCD results |J. These two approaches yield remarkably consistent 
determinations of \V u b\, but have large theoretical uncertainties.^ 



#Talk presented at the XXXIIIrd Rencontres de Moriond: Electroweak Interactions and Unified 
Theories, Les Arcs, Savoie, France, March 14-21 1998. 

1 V u b is one of the least precisely known parameters of the standard model. The other poorly known 
CKM matrix element, Vtd, is related to V u b through the unitarity triangle. 
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To reduce these uncertainties, it will be advantageous to consider different observables. 
The problem is that there is almost no overlap between quantities sensitive to \V u b\ which can 
be reliably calculated theoretically and those which can be measured experimentally. Here I 
discuss two proposals which bridge this gap to some extent. In Section II the possibility of 
extracting \V u b\ from the hadron invariant mass spectrum in inclusive semileptonic B decay is 
reviewed. Section III concerns extracting \V u b\ from a double ratio of form factors in exclusive 
semileptonic B and D decays to p and K*. Section IV contains some conclusions. 



II. V ub FROM INCLUSIVE B DECAYS 



The traditional method for extracting \V u b\ from experimental data involves a study of the 
electron energy spectrum in inclusive charmless semileptonic B decay @. In the B rest frame, 
electrons with energies in the endpoint region E e > {m 2 B — m 2 D )/2mB must arise from b —>■ u 
transitions. There has been considerable theoretical progress recently in our understanding of 
inclusive semileptonic B decay @-§|, based on the use of the operator product expansion (OPE) 
and heavy quark effective theory. At leading order in the AqcdA 7 ^ expansion the B meson 
decay rate is equal to the b quark decay rate. There are no nonperturbative corrections of order 
Aqcd/^6- In the electron endpoint region our calculational ability is lost, since the size of this 
region m 2 D /2mB — 330 MeV is comparable to (rriB — rrib)/2. An infinite set of higher order 
terms in the OPE, which extend the endpoint from mb/2 to wb/2, yield singular contributions 
to dr/d-E e that are equally important integrated over such a small endpoint region. 

In the future, it may be possible to determine \V u b\ from a comparison of the measured 
hadronic invariant mass spectrum in the region sh < rn 2 D with theoretical predictions 
Here Sh = {pb — q) 2 , where Pb is the B meson four- momentum, and q = p e + pp is the sum 
of the lepton four-momenta. An obvious advantage to studying this quantity rather than the 
lepton energy spectrum is that most of the B —>■ X u e v decays are expected to lie in the region 
Sh < raj), while only a small fraction of the B — > X u ev decays have electron energies in the 
endpoint region. Both the invariant mass region, sh < mfj, and the electron endpoint region, 
E e > (m 2 B — m 2 D )/2mB, receive contributions from hadronic final states with invariant masses 
that range up to m^. However, for the electron endpoint region the contribution of the states 
with masses nearer to is kinematically suppressed. This region is dominated by the ir and 
the p in the ISGW model ||10|| , with higher mass states making only a small contribution. The 
situation is very different for the low invariant mass region, sh < rn 2 D . Now all states with 
invariant masses up to contribute without any preferential weighting towards the lowest 
mass ones. In the ISGW model the 7r and the p mesons comprise only about a quarter of the B 
decays to states with sh < fnjj- Consequently, it is much more likely that the first few terms in 
the OPE will provide an accurate description of B semileptonic decay in the region sh < tt? d 
than in the endpoint region of the electron energy spectrum. In fact, from a theoretical point 
of view, the cut sh < fn 2 D provides the optimal kinematical separation between inclusive b — ► u 
and b — > c decays. A modest cut on the electron energy, which will probably be required 
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experimentally for the direct measurement of sh via the neutrino reconstruction technique, 
will not destroy this conclusion. 

To begin with, consider the contribution of dimension three operators in the OPE to the 
hadron mass squared spectrum in B — ► X u e v decay. This is equivalent to b quark decay and 
implies a result for dT/dE ds (where E = p b ■ (p b — q)/m b and s Q = (p b — q) 2 are the energy 
and invariant mass of the strongly interacting partons arising from the b quark decay) that can 
easily be calculated using perturbative QCD up to order a 2 (3o- Even at this leading order in 
the OPE there are important nonperturbative effects that come from the relation between the 
b quark mass and the B meson mass, m B = m b + A + 0(Aq CB / m b ). The most significant effect 
comes from A, and it relates the hadronic invariant mass sh to Sq and E via 

s H = s + 2AE + A 2 . (1) 

Changing variables from (s , E ) to (s#, E ) and integrating E over the range 

y/s^ — A < E < — (sh — 2Am B + m 2 B ), (2) 

gives dT/dsn, where A 2 < sh < rn B . Feynman diagrams with only a u-quark in the final state 
contribute at s = 0, which corresponds to the region A 2 < s# < Am#. 

Although dT/dsn is integrable in perturbation theory, powers of a s \n 2 [(sH — Aiiie) /rn 2 B ] 
occur in the invariant mass spectrum. This shows that perturbative and nonperturbative 
corrections are both important for sjj < Am^. (In the m b — > oo limit perturbative corrections 
are important in a slightly larger region since a s ln(%/m|) ~ 1 for sh ~ Amg.) While dT/dsn 
cannot be reliably predicted for sh ~ Atob, the behavior of the spectrum for sh ~ Arris 
becomes less important for observables that average over larger regions of the spectrum, such 
as dT/dsn integrated over sh < A 2 , with A 2 significantly greater than Arris- 

In Fig. 1 we plot the quantity f (A 2 , A) defined by 

a 2 , dT(B^X u eu) G 2 F m% lTrl2 f A\ 5 ~ 2 - N 

ds H^~ A — — = 1§T? Kfc 2 i r(A , A) , 3 

o ds H 1927r d V m B J 

as a function of A 2 for A = 0.2, 0.4 and 0.6 GeV in the range Am# < A 2 < 4.5 GeV 2 , including 
terms up to order a 2 /3 (using a s {m b ) = 0.2). These curves approach f(m 2 3 ,A) ~ 0.73 as 
A 2 — ► m 2 B \\11\\ . f (A 2 , A)/f (m 2 B , A) is the fraction of events with hadronic invariant mass less 
than A 2 . It is mostly the ability to compute f(A 2 ,A) and our knowledge of the value of A 
which determine the uncertainty from theory in a value of \V ub \ extracted from the invariant 
mass spectrum in the region sh < A 2 . 

In the low mass region, sh ~ Am#, nonperturbative corrections from higher dimension 
operators in the OPE are very important. Just as in the case of the electron spectrum in 



the endpoint region [12], the most singular terms can be identified and summed into a shape 



function, S(sh)- Neglecting perturbative QCD corrections, we write 
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1.5 2 2.5 3 3.5 4 4.5 
A 2 [GeV 2 ] 

FIG. 1. The function f(A 2 ,A) defined in Eq. (|) as a function of A 2 for A = 0.2 GeV (dotted 
curve), 0.4 GeV (solid curve), and 0.6 GeV (dashed curve). Note that m 2 D = 3.5 GeV 2 . 

It is convenient to introduce the scaled variable y = %/Amf,. Then |7j 

s(y) = E MtP h 2 v n+2 ( 3 - 2 v) *a - y)] > ( 5 ) 



n=0 



n\ A" dy n 



where S(y) = Arrib S(sh) is dimensionless. The matrix elements A n are the same ones that 
determine the shape functions for the semileptonic B decay electron energy spectrum in the 
endpoint region, and also the photon energy endpoint region in weak radiative B decay, 

(B(v)\ iD^ . . . iD^ n \B(v))/2m B = A n v lll ... + (terms involving g^.) . (6) 

The AnS have dimension of [mass] n , and hence the coefficients A n /A n are dimensionless num- 
bers of order one. The first few A n 's are A — 1, A\ — 0, A 2 = — Ai/3, etc. 

The shape function S(y) is an infinite sum of singular terms which gives an invariant mass 
spectrum that leaks out beyond y = 1 (i.e., sh = Arrib)- For y ~ 1, all terms in Eq. (^) are 
formally of equal importance. Since Am?, ~ 2 GeV 2 is not too far from m 2 ^, it is necessary to 
estimate the influence of the nonperturbative effects on the fraction of B decays with invariant 
hadronic mass squared less than A 2 . It is difficult to obtain a model-independent estimate of 
the leakage of events above an experimental cutoff sh = A 2 , given that we can estimate only 



the first few moments, A n . However, in the ACCMM model JT3| with reasonable parameters, 
the shape function S(y) causes a small (i.e., ~ 4% with A = 0.4 GeV, and perturbative QCD 
corrections neglected) fraction of the events to have sh > rn 2 D 0. Moreover, this leakage 
depends primarily on A, and only to a lesser extent on other ingredients of the model. 

Thus the analysis of both perturbative and nonperturbative corrections implies that the 
uncertainty in the determination of \V u b\ from the hadronic invariant mass spectrum in the 
region sh < itl 2 d is largely controlled by the uncertainty in A, or equivalently, by that in the b 
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quark mass. (There is a so-called renormalon ambiguity in A and in the b quark pole mass. For 
the physically measurable quantity in Eq. (^), this is cancelled by a similar ambiguity in the 
perturbative series in T.) To measure \ V u b\ with a theoretical uncertainty below ~ 10%, A has 
to be determined [0 (with better than ~ 100 MeV uncertainty), and the cut on sh has to be 
as close to m 2 D as possible. If the experimental resolution forces one to consider a significantly 
smaller region then the theoretical uncertainties will be larger. 



III. V ub FROM EXCLUSIVE B DECAYS 

Heavy quark symmetry [15|] is much less predictive for heavy to light decays than it is for 



heavy to heavy transitions. In the infinite mass limit not all form factors are related to one 
another, and their normalization is not fixed at any kinematic point. There are still relations 



between semileptonic B and D decays to the same charmless exclusive final state |16| , such as 
between B — > p£v and D — > p£u, or between B — > tt£u and D — > tt£v. The order l/m c j, 
corrections to the infinite mass limit may be sizable, however, and for final state pions there are 
additional complications since m w is comparable to the mass difference between the vector and 



pseudoscalar mesons ||17|| . The question is whether we can construct an observable sensitive to 



V u b that is (almost) free of \jm c h corrections. 

The basic idea [0,0 is to use heavy quark symmetry to relate the SU(3) violation between 
D — > K* £ v and the Cabibbo suppressed decay D — > p £ v to those that occur in a comparison 
of B — > K*£ I (or B — > K* v v) with B — > pi v. Then experimental data on B — ► K*£ £ 
in conjunction with data on D — > piv and D — > K* Iv can be used to determine \V u b\. 
This proposal is complementary to other approaches for determining \ V u b\, since it relies on the 
standard model correctly describing the rare flavor changing neutral current process B — ► K *£ £. 

We denote by g^ H ^ v \ f ^ H ~~* V \ and a±^ the form factors relevant for semileptonic tran- 
sitions between a pseudoscalar meson containing a heavy quark H (H = B,D), and a member 
of the lowest lying multiplet of vector mesons V (V — p, K*, lu), 



(V(p', e) I qj, Q \H(p)) = i g^ V) e, uXa f (p + p') x (p - V 'Y , (7) 
(V(p', e) | q 7,75 Q \H(p)) = f {H ^ v) e£ + af ~> V) (e* • p) (p + p')» + a {H ^ V) (e* ■ p) (p - p'), . 



We view the form factors as functions of the dimensionless variable y — v ■ v', where p = m# v, 
p' = my v', and q 2 = (p — p') 2 = m 2 H + m 2 / — 2m# my y. (Although we are using the variable 
v ■ v', we are not treating the quarks in V as heavy.) Assuming nearest pole dominance for the 
q 2 dependences, the D — > K* £ v form factors are ||T9| 



f (D^K-),, A = (1-9 ± 0.1) GeV 
W) 1 + 0.63 {y - 1) 

u/) = 



,(D-*K*) fiii>l _ (0-18 ± 0.03) GeV" 1 

1 + 0.63 (y-l) ' 
(^) M = (0-49 ± 0-04) GeV' 1 

^ iyj 1 + 0.96 (y-l) ' lJ 
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The shapes of these form factors are beginning to be probed experimentally ||19|| . The form 
factor a_ is not measured because its contribution to the D — > K* t v decay amplitude is 
suppressed by the lepton mass. These form factors are measured over the kinematic region 
1 < y < {m 2 D + m 2 K *) / (2rriD m K *) ~ 1.3. Note that f(y) changes by less than 20% over this 
range. The full kinematic region for B — > ptv is much larger, 1 < y < 3.5. In the following 
analysis we will extrapolate the measured D — > K* form factors to 1 < y < 1.5. (The validity 
of this extrapolation can be tested . ) 

The differential decay rate for semileptonic B decay (neglecting the lepton mass, and not 
summing over the lepton type t) is 



dT(B -f ptv) 



r<2 



V, 



lib 



dy 

Here S( H ^ v \y) is the function 



48 7T 3 



m B 



mlS^Ky) 



(9) 



y 2 -i 

+ 4Re 



f^ v \y) (2+y 2 -6yr + 3r 2 ) 
a^ V \y) f {H - V) (y)]m 2 H r (y - r){y 2 - 1) 



+ 4 



(y) m\ir 2 {tf 



1) 2 + 



(y) 



m H r 



'l + r 2 -2yr)(y 2 -l) 



(10) 



y 2 - 1 f {H ^ v \y) (2 + y 2 - 6yr + 3r 2 ) [1 + 6^ v \y)] , 

with r = my /ran- The function $( H ~* V '> depends on the ratios of form factors 
and g( H ~^ v ) / f( H ~^ v \ S( B ~* p \y) can be estimated using combinations of SU(3) flavor symmetry 
and heavy quark symmetry. SU(3) symmetry implies that the B° — > p + form factors are equal 
to the B — > K* form factors and the B" — > p° form factors are equal to l/y/2 times the B — > i^* 
form factors. Heavy quark symmetry implies the relations [|16|] 



(/, a + , 9 r-^ 



fm B s 


r 


~a s {m h y 


-6/25 






_a s (m c )_ 





(/, a+, 9) 



(11) 



where we used a^~* K ^ = —a^^"^ ' , valid in the large m c limit. 

Using Eq. flTT| ) and S77(3) symmetry to get B° — > p + tv form factors (in the region 1 < y < 
1.5, corresponding to q 2 > 16 GeV 2 ) from those for D — > K*tv given in Eq. (§) yields S (B ^ p) (y) 
plotted in Fig. 1 of Ref. Hi~8~| |. (The numerical values in Eq. (§) differ slightly from those used 
in Ref. ||18[| .) This prediction for S^ B ^ 9 ^ can be used to determine \V u b\ from the B ^ ptv 
semileptonic decay rate in the region 1 < y < 1.5. We find that about 20% of B° — > p + tv 

This method is 
Typical 



(D—>K*) 



y<1.5 



5.9 IK 



ub 



decays are in the range 1 < y < 1.5, and B(B° — > p + £z/ 

model independent, but cannot be expected to yield a very accurate value of \V u b 
SU (3) violations are at the 10 — 20% level and similar violations of heavy quark symmetry are 
expected. In this region \^ B ~* p \y)\ defined in Eq. ( |10"D is less than 0.06, indicating that a^ 3 ^^ 
and g( B ~* pS > make only a small contribution to the differential rate in this region. Thus the main 
uncertainties are SU{3) and heavy quark symmetry violations in the f( H ~* v *> form factor only; 
these are precisely the ones we can eliminate. 
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FIG. 2. Feynman diagram that gives the leading contribution to R(l) — 1. The dashed line is a ir 
or an r\. The black square indicates insertion of the weak current. 

Ref. proposed a method for getting a value of S( B ~* p '(y) with small theoretical uncer- 



tainty using a "Grinstein-type" ETJ double ratio 



R(y) = [f iB ^ p) (y)/f (B ^ ) (y)]/[f (D -* p) (y)/f iD ^ ) (y)} , (12) 

which is unity in the limit of SU(3) symmetry or in the limit of heavy quark symmetry. 
Corrections to the prediction R(y) = 1 are suppressed by m s /m c ^ (m Uj d <C m s ) instead of 
7?t, s /Aqcd or Aqcd/^c,6- The leading deviation of -R(l) from unity arising from the Feynman 
diagrams in Fig. 2 has been estimated using chiral perturbation theory. These yield a calculable 
non-analytic ^m q dependence on the light quark masses, and such corrections cannot arise from 
other sources. The result is R(l) = 1 — 0.035 (7 (72 j22fl , where g is the DD*n coupling and g% is 
the pun coupling. Experimental data on r — > u tt v t decay gives #2 — 0.6 ||23|| . Estimates of 
g vary between near unity and much smaller values ||24j| . There may be significant corrections 



to R(l) from higher orders in chiral perturbation theory. However, the smallness of our result 
lends support to the expectation that -R(l) is very close to unity. There is no reason to expect 
any different conclusion over the kinematic range 1 < y < 1.5. 
Since R(y) is very close to unity, the relation 

s^'(y) = s^«-\y) 2 fj^yy _ (13) 

/ \V) \m B -m K *J 

and measurements of \f( D ~* K ** > \, \f^ D ^ p ^\, and S^ B ^ K *^ will determine S^ B ^ p,) with small theo- 
retical uncertainty. The last term on the right hand side makes Eq. ( p.3[ ) equivalent to Eq. (|T2]) 
in the y — > 1 limit. The ratio of the (2 + y 2 — 6yr + 3r 2 ) [1 + 5^ B ^ v \y)] terms makes only 
a small and almost y-independent contribution to S^ B ^ P ^ / S^ B ^ K *^ in the range 1 < y < 1.5. 
Therefore, corrections to Eq. (|H|) are at most a few percent larger than those to R(y) = 1. 

\f^ K *~>\ has already been determined. |/ < - D ^ p - ) | may be obtainable in the future, for exam- 
ple from experiments at B factories, where improvements in particle identification help reduce 
the background from the Cabibbo allowed decay. The measurement B(D — > p° £v)jB(D — > 
K*°£u) = 0.047 ± 0.013 |5[ already suggests that \f( D ^p) / f( D ^ K *)\ is close to unity. As- 
suming 577(3) symmetry for the form factors, but keeping the explicit my-dependence in 
S^ D ~^ v '(y) and in the limits of the y integration, the measured form factors in Eq. (|8|) im- 
ply B(D -> p°£u)/B(D -> K*°£u) = 0.044 @. 

cj(B—tK*) j g obtainable from experimental data on B — > K*££ or B — > K* vv. While the 
latter process is very clean theoretically, it is very difficult experimentally. A more realistic goal 
is to use B — > K*££, since CDF expects to observe 400 — 1100 events in the Tevatron Run II (if 
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the branching ratio is in the standard model range) ||26|| . The uncertainties associated with long 



distance nonperturbative strong interaction physics in this extraction of S^ B ~* K *'(y), averaged 
over the region 1 < y < 1.5, are probably less than 10% [p2| . Consequently, a determination 
of \V u b\ from experimental data on D — > K*£u, D —>■ p£v, B — ► K*££ and B ^ p£u with an 
uncertainty from theory of about 10% is feasible. If a precise value of \V u b\ is available before 
B — > K*££ is measured, then we get an accurate standard model prediction for the B — > K*££ 
decay rate in the region 1 < y < 1.5. Comparison with data may signal new physics or provide 
stringent constraints on extensions of the standard model. 

IV. CONCLUSIONS 

The present determinations of \V u b\ rely on comparing experimental data with model cal- 
culations, and therefore suffer from theoretical uncertainties of order 30%. (This is hard to 
quantify, and such a number is necessarily ad hoc.) To reduce these uncertainties one needs 
to consider somewhat different observables for which the theoretical predictions are less model 
dependent than those for the endpoint region of the inclusive electron spectrum and for the 
total exclusive B — > 7r Iv or B — ► p£v decay rates. In this talk I reviewed two ideas which 
seem promising to me: i) extracting \V u b\ from the hadronic invariant mass spectrum in inclu- 
sive semileptonic B decays; and ii) using heavy quark and chiral symmetries for form factors of 
exclusive semileptonic B and D decays to vector mesons. These may lead to model independent 
determinations of \V u b\ with an uncertainty from theory of about 10%. Lattice calculations p7| 
and dispersion relation constraints on form factors |2(J will also be important. 

There is not really one "gold-plated" observable for extracting \V u b\. To reduce the strong 
interaction model dependence, several measurements will be needed to guide us which approx- 
imations and expansions have smaller uncertainties. At the 10% level consistency between 
different determinations of \V u b\ will be necessary to have confidence that the uncertainties are 
indeed so small. I am hopeful that this will be achieved within the next few years. 
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